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Abstract 

We study a single statistical amphiphilic copolymer chain AB in a selective solvent 
(e.g. water). Two situations are considered. In the annealed case, hydrophilic (A) and 
hydrophobic (B) monomers are at local chemical equilibrium and both the fraction of A 
monomers and their location along the chain can vary, whereas in the quenched case (which 
is relevant to proteins), the chemical sequence along the chain is fixed by synthesis. In both 
cases, the physical behaviour depends on the average hydrophobicity of the polymer chain. 
For a strongly hydrophobic chain (large fraction of B), we find an ordinary continuous 
9 collapse, with a large conformational entropy in the collapsed phase. For a weakly 
hydrophobic, or a hydrophilic chain, there is an unusual first-order collapse transition. In 
particular, for the case of Gaussian disorder, this discontinuous transition is driven by a 
change of sign of the third virial coefficient. The entropy of this collapsed phase is strongly 
reduced with respect to the 9 collapsed phase. 
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Copolymeres aleatoires hydrophiles-hydrophobes 

Resume 

Nous etudions un copolymere aleatoire amphiphile AB dans un solvant selectif (par 
exemple, de I'eau). Nous considerons deux cas. Dans le cas du desordre mobile, les 
monomeres hydrophiles (A) et hydrophobes (B) sont a I'equihbre chimique local, et la 
fraction de monomeres A ainsi que leur position dans I'espace peuvent varier, alors que 
dans le cas du desordre gele (qui est relie au probleme des proteines) , la sequence chimique 
est fixee par synthese. Dans les deux cas, le comportement de la chaine depend de son 
hydrophobicite moyenne. Pour une chaine fortement hydrophobe (grande fraction de B), 
on trouve un point d'effondrement 6 continu ordinaire, avec une grande entropie conforma- 
tionnelle. Pour une chaine faiblement hydrophobe ou hydrophile, on trouve une transition 
inhabituelle du premier ordre. En particuler, dans le cas du desordre gaussien, cette tran- 
sition discontinue est pilotee par un changement de signe du troisieme coefficient du viriel. 
L'entropie de cette phase coUapsee est fortement reduite par rapport a celle d'un point 9 
ordinaire. 
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I. Introduction 

The study of random heteropolymers in solutions has been recently developped, the 
main reason being a plausible connection of this problem with the protein folding puzzle 
Various kinds of quenched randomness have been considered. Here, we wish to 
emphasize the role of the solvent (water) in the folding process. Indeed, it is widely 
believed that the hydrophobic effect is the main driving force for the folding transition: 
in a native protein, the hydrophobic residues are hurried in an inner core, surrounded by 
hydrophilic (or less hydrophobic) residues. In this paper, we consider a simple model, 
where the monomers of a single chain are (quenched) randomly hydrophilic or hydrophobic. 
They interact with the solvent molecules through an effective short range interaction, 
which is attractive or repulsive, depending on their nature. A related model has been 
considered by Obukhov but he did not study its phase diagram. We first consider 
the corresponding annealed case, which is simpler and gives useful insights into the nature 
of the low temperature phase. It may well be of interest for hydrophilic polymers in 
the presence of complexing amphiphilic molecules. In section II, we describe the model. 
Section III deals with the annealed case, and section IV with the quenched case, both for 
Gaussian and binary distributions of hydrophilicities. 

II. The model 

We consider a chain of N monomers which can be either hydrophilic (L) or hydropho- 
bic (P). The hydrophilicity degree of monomer i is measured by A^, and the monomer- 
solvent interaction is assumed to be short-ranged ((^-interaction). More precisely, we con- 
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sider a monomer-solvent interaction of the form: 

N+l Af 
i=l a=l 

where the fi are the positions of the N monomers and the Ra are the positions of the J\f 
solvent molecules ^. Here, the Aj are independent annealed or quenched random variables. 
A positive ( resp. negative) A corresponds to a hydrophilic (resp. hydrophobic) monomer. 
Using the incompressibility condition of the chain-solvent system i.e. the sum of the 
monomer and solvent concentration is constant: Pm(r) + Ps(r) — POi solvent degrees 
of freedom can be eliminated, and the full Hamiltonian reads : 

+ ^ Yl '^sSifi- Tj) 5{fj - fk) (2) 

+ ^ 5Z S{fi - fj) S{fj - fk) S{fk - n) 

In equation (2), the vq term accounts for entropic excluded volume interactions as well as 
effective AB interactions described by Flory parameter xab ^- Usually, when considering 
a polymer collapse transition, it is sufficient to include terms up to w^. Here, we have 
included a repulsive four-body term W4 which will play an essential role. 
In the continuous limit the partition function reads: 

ZiXis)) = J Pf(.)exp (-^J\s (^^y - pniXis))\ (3a) 



^ Note that —X/kT is identical to the standard Flory parameter x, often used in the 
polymer literature. 

^ In the Flory-Huggins model of polymer solutions, vq — a'^{l — 2xab), ws = and 
W4^ = a^, where a is the monomer size^^^. 



with 

PniX{s)) = - / ds ds' {vo + /3(A(s) + X{s'))) 5{f{s) - f{s')) 
^ Jo Jo 

pN pN pN 

+ ^ / ds ds' ds" 5{f{s)-f{s')) 5{f{s')-f{s")) 
6 Jo Jo Jo 

pN pN pN pN 

+ S / ^^"/ (^s'''<^(r(s)-f(s0)<^(^(s0-f(O)<^(^(s'0-^(O) (3&) 

24 Jo Jo Jo Jo 

where we have imphcitly assumed that the origin of the chain is pinned at and its 
extremity is free. The parameter a denotes the statistical segment (Kuhn) length , and 
d denotes the dimension of space. In the following, we shall study the case where the 
hydrophilicities X{s) are random independent Gaussian variables, with average Aq and 
variance A: 

We shall also discuss briefly the technically more complicated case where the Aj are random 
independent binary variables, taking the value A^ with probability p for a hydrophilic 
monomer and Xb with probability q = 1 — p for a hydrophobic monomer. This case is 
more relevant to synthetic random copolymers. 

The annealed case amounts to take the average of the partition function Z on the 
randomness, whereas the quenched case corresponds to taking the average of log Z. Since 
the quenched case is achieved by introducing replicas and using the identity logZ = 
lim„_»o ~^ , we compute directly the average for integer n, and then treat separately 
the annealed (n = 1) and quenched (n = 0) cases. For Gaussian disorder and integer n, 
the average yields: 
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exp 

with 



(o2\2 pN pN pN n \ 

^ dsj^ ds' ds" J] 5(f«(s)-r;(/))5(f,(s)-f,(0)l (5a) 



,{s)-fa{s')) 
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Ag = T^ivo + 2/3 Ao) / ds / ds' V5(ra( 

„Ar „iv 

/ ds ds' ds" y25{fa{s)-fa{s'))5{fa{s')-fa{s")) 

Jo Jo Jo 



fN pN pN pN 

+ / ds ds' / ds" / ds'" 



Jo Jo Jo 



24 

n 

Yl ^i^ais) - fa{s')) S{fa{s') - ra{s")) SiVais") - fa{s"')) (56) 



a=l 

whereas for binary disorder: 



with 



Ab = ^ ds ds'V] Sifais) - fais')) 

^ a=l 

+ / '^"^ / ds' I ds" ^(ra(g) - ra{s')) 5{fa{s') - fa{s")) 
" Jo Jo Jo 

pN pN pN pN 

+ ^ / ds ds' / ds" / ds'" 



JO JO JO 



24 

n 

E '^(^"'«(^) - ^ais')) Sifais') - fais")) Sifais") - fais'")) 



a=l 



(66) 



We now study separately the annealed and quenched cases. 
III. The annealed case 

We set n = 1 in equations (5) and (6). Denoting the monomer concentration by p(r), 
we assume the validity of ground state dominance ^^'^^ . The ground state wave-function 



(fi{r) is normalized, and we perform a saddle-point approximation of the partition function. 
We quote here only the results, leaving the details of the calculations for the quenched case 
(section IV). These approximations yield the following annealed free energy per monomer: 
a) Gaussian disorder 

+N^'^ J d'^r ip%f) + J d'^r ip\r) - Eq J d'^r (p{r)'' (7a) 

where 

w'^ = W3- SP^X^ {7b) 

The energy Eq appearing in (7a) is the Lagrange multiplier enforcing the normalisation of 
the wave-function (p. The equation of state is obtained by minimizing the free energy (7a) 
with respect to the wave-function (f{r). The monomer concentration is given by: 

p{f^ = N<f\r) (8) 

The phase diagram of the system can be easily discussed: for any non-zero concen- 
tration of hydrophobic monomers, there is a transition between a swollen phase (p = 0) at 
high temperature, and a collapsed phase (finite p) at low temperature. The nature of the 
transition depends on the average hydrophobicity of the chain Aq and variance A. 

For strong average hydrophobicity (Aq < 0), there is a second order collapse transition, 
similar to an ordinary 9 point for polymers in a bad solvent ^^^^ . This occurs when the 
coefficient of the cp'^ term in (7a) becomes negative while the coefficient of the cp^ term is 



still positive. This occurs if the average hydrophobicity Aq satisfies the condition: 



When condition (9a) is not satisfied, the collapse transition becomes first order, and is 
driven by the change of sign of w'^, i.e. by a change of sign of the effective third virial 
coefficient. This is quite an unusual mechanism. 

Such is the case if the system is weakly hydrophobic, or hydrophilic. The minimization 
of the free energy (7a) leads to a non-linear partial differential equation for (p which cannot 
be solved analytically. It is possible to restrict the minimization to a subset of Gaussian 
wave- functions, and this will be discussed in the section on quenched systems (see eq.(20) 
below). Within this Gaussian approximation, the first order transition temperature T* is 
given by the equation: 




(9a) 



and the transition temperature is given by: 



2Ao 



(96) 




(9c) 



b) Binary disorder 



In this case, the free energy per monomer reads: 




■) 



(10a) 



where 
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By taking derivatives with respect to the chemical potentials A±, it is easily seen that the 
concentrations of hydrophilic and hydrophobic monomers are given by: 



p+(r-) = iV(^2(r) 



p_(r) = 7V^^(r) 



(11) 



Note the Boltzman weights in (11), which govern the chemical composition of the polymer. 
It clearly appears that the regions of high monomer concentration are hydrophobic, whereas 
those of low concentration are hydrophilic. The main conclusions of the Gaussian case 
above remain valid, except that the transition temperatures cannot be calculated explicitly. 
We find again an ordinary second order 9 point for strongly hydrophobic chains, and a 
first-order collapse for weakly hydrophobic or hydrophilic chains. 

As can be seen by a numerical study of this case, the transition is strongly first-order, 
and this implies the existence of a large latent heat at the transition, and thus a sharp drop 
of the entropy. This drop in entropy can be interpreted as the formation of a hydrophobic 
core, surrounded by a hydrophilic coating, which is energetically favourable, but entropi- 
cally unfavourable. This is confirmed by equations (11), which show that high concentra- 
tion regions are hydrophobic, whereas low concentration regions are hydrophilic. Finally, 
we see from (11) that at low temperature, the total number of hydrophilic monomers goes 
to zero, whereas the number of hydrophobic monomers goes to N. 

IV. The quenched case 

We have to take the n = limit in (5) and (6). Introducing the parameters qab{r,f') 
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, with a < b, and Pa(r) by: 

Qahif, r)= ds S{fa{s) - r) 5{n{s) - f 
Jo 



.N (12) 

Pair) = ds S{ra{s) - r) 
Jo 



we may write equations (5) as: 
a) Gaussian disorder 



Zg = J T^Qabir, f')Vqabir, r')Vpa{r)V(j)a{r) exp {G{qab-, Qab, Pa, (i>a) + logC(9a6, (i>a)) 

(13a) 

where qabif^f) and (/>a(r) are the Lagrange multiphers associated with (12), and: 

G{qab, qab, Pa, </>a) = / E i^'Pai^^^ " (^0 + W^o) ^ - ^P^(0 - ^P^(r)) 

+ dr dr'^{iqab{r,r')qab{r,r') + qab{r,r')Pa{r)Pb{r')) 

a<b 

(136) 



and 



a 

xexp^-^y ^ 1^^^ -ij ds ^(j)a{ra{s)) -i ds^^qab{ra{s),h{s)) 

(13c) 

where w'^^ is defined in equation (7b). 

Using the identity of equation (13c) with a Feynman integral , we can write: 

Q{qab,(t^a)^ I n^^'^a<ri...r;|e-^^"(^<''"'^'^)|0...0> (14a) 

a 
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where is a "quantum-like" n — > Hamiltonian, given by: 

= + J]^</'a(ra) + J2iqab{fa,fb) (146) 

a a a<b 

Anticipating some kind of (hydrophobically-driven) collapse, we assume that we can use 
ground-state dominance to evaluate (14), and write, omitting some non-extensive prefac- 
tors: 

/ \ (15) 

= exp -N min {< *|i^n|* > -Eq (< > -1)} 

where Eq is the ground state energy of H^. At this point, the problem is still untractable, 

and we make the extra approximation of saddle-point method (SPM). The extremization 

with respect to qab reads: 

iqabir. n = -P^ X^Paif) Pbin (16) 

This equation shows that replica symmetry is not broken, since qab is a product of two 
single-replica quantities pa- To get more analytic information, we use the Rayleigh-Ritz 
variational principle to evaluate Eq. Due to the absence of replica symmetry breaking 
(RSB), we further restrict the variational wave- function space to Hartree-like replica- 
symmetric wave-functions, and write: 

n 

^{fi,...,fn)^l[^{fa) (17) 
a=l 

Because of replica symmetry, we can omit replica indices and take easily the n ^ limit. 
The variational free energy now reads: 

-PFhiq, g, p, 0, ^) = y dr (^ip(r)</>(r) - (^o + 2/3Ao) ^ - "fp^f) - ||p'(r)) 
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— - J '^^ J dr' (iq{f,f^)q{f,r) + P'^X^ q{f,r)p{f^p{r)^ 

+NEo Q dV'(^) - (18) 

The SPM equations read: 
p(r) = Nip^{r) 

q{r,r') ^ Nip^{r)<p^{r') 

r (19) 

i<Pif) = {vo + 2/3Ao)p(r) + ^p\f) + ^p^f) + (^^>? J d''r'q{f, r')p{r') 

iq{r, r') = -I3'^\^p{7^p{r') 
We still have to minimize with respect to the normalized wave-function ip{f). This leads 
to a very complicated non-linear Schrodinger equation, and we shall restrict ourselves to 
a one-parameter family of Gaussian wave-functions of the form: 



(20) 



where R is the only variational parameter. 

Using equations (18) and (19), the variational free energy per monomer (/g = Fq/N) 
reads: 

-Eo J d'^r^^r) (21) 

Note that the only difference between this free energy and the annealed one of eq.(7a), is 
the presence, in the quenched case, of an additional disorder induced term, which we have 
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grouped together with the term for reasons that will become clear below. Using the 
Gaussian wave function (20), and the value of w'^ given in (7b), the free energy reads: 

r,-r- 1 {vo + 2/3 Ao) N 



+ 



8i?2 (20F)^ 2 Rd 

'—^(3-^/^-2-^)] i^f+f^Y'' ^ f^Y (22) 
Try 2 ^ V ^R'^ V(327r3); 24 Uv ^ ^ 



X^ttVsY 6 (27 

At low temperatures, one has to study the sign of the third term of eq(22) . This sign 
indeed becomes negative at low temperature, due to disorder fluctuations, yielding the 
following results: 

i) Ao > : the hydrophilic case. 

In our approximation, we find that there is a first order transition towards a collapsed 
phase (i? ~ N^/'^) induced by a negative three-body term (and stabilized by a positive 
four-body term). This transition is neither an ordinary 9 point (since the two-body term 
is positive), nor a freezing point (the replica-symmetry is not broken). Note that since our 
approach is variational, the true free energy of the system is lower than the variational one, 
and we thus expect the real transition to occur at an even higher temperature. Since the 
transition is first order, we expect metastability and retardation effects to be important 
near the transition; we also expect that (due to the latent heat), there will be a reduction 
of entropy in the low temperature phase, compared to an ordinary second order 6 point. 

ii) Ao < : the hydrophobic case. 

In this case, defining the two temperatures: 

To = 2|Ao|/^;o 

T, = (23) 
7103(1 -3^^/2/2'^) 
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leads to two possible scenarii, in a very similar way as in the annealed case: 

a) Weakly hydrophobic case : Ti > Tq 

The collapse transition is again driven by the disorder fluctuations of the three-body 
interactions. The resulting first-order transition is very similar to the hydrophilic case i). 

b) Strongly hydrophobic case : Tq > Ti 

The collapse transition is now driven by the strong two-body Aq term. The resulting 
phase transition is very similar to an ordinary 9 point, and is therefore second-order. 

The phase diagram of the quenched case seems quite similar to the annealed case, 
but, as we shall see on the case of binary disorder, the physics is quite different. Note also 
that the transition temperature in the quenched case is lower than in the corresponding 
annealed case, as can be checked from (7a) and (21). This may be viewed as an effect of 
geometrical frustration induced by the chain constraint, 
b) Binary disorder 

Using the same approximations as in the previous case, namely ground state domi- 
nance, no replica symmetry breaking and Hartree variational wave- function, the free energy 
per monomer reads: 

+NPXa J d''np\r) - ^^^7— (^)' log (/ dWif^e^'f^^' ^'^^^^^ (24) 



where A = -|- A^. 

In this case, we see that even the Gaussian variational form is not tractable, but a 
numerical study shows that the system has exactly the same phase diagram as for Gaussian 
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disorder. However, the interest of this model is that it allows to calculate the concentration 
of hydrophilic and hydrophobic monomers. We obtain: 



and thus (as expected for a quenched chain) , the overall chemical composition of the chain 
remains unchanged. 
V. Conclusion 

We have studied a simple model of a randomly hydrophilic-hydrophobic chain in 
water. The randomness may be annealed or quenched. The low temperature phase is 
always collapsed; depending on the degree of hydrophobicity of the chain, this collapsed 
phase is reached either through an usual second order 9 point, or through a first order 
transition. In all cases (see eq.(ll) and (25)), the interior of the compact globule is mainly 
hydrophobic, whereas its exterior is hydrophilic (as expected). In the case of a first order 
transition, the existence of a latent heat implies a decrease of the entropy of the collapsed 
chain, due to the burial of the hydrophobic monomers in the core. Studies of similar 
models for semi-dilute solutions and melts in random or block copolymers are currently 
under way. 




(25) 



It is easily seen on (25) that the two concentrations satisfy the sum rules: 




(26) 
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